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Layers  in  the  Atmosphere 

ABSTRACT 

Measurements  were  made  of  the  scattering  properties  of  thin  turbulent 
layers  at  and  above  the  tropopause.  The  Millstone  Hill  L-band  radar  was 
used  to  measure  the  backscatter  cross  section  per  unit  volume  of  these  layers 
as  a  function  of  time  and  space.  An  X-band  forward  scatter  link  was  set  up 
between  Wallops  Island,  Virginia  and  Westford,  Massachusetts  to  observe 
scattering  from  these  layers.  Although  the  radar  could  not  provide  observa¬ 
tions  of  the  common  volume  of  the  forward  scatter  link,  for  days  where  no 
clouds  were  observed  in  the  vicinity  of  the  tropopause,  the  radar  observations 
of  layers  near  the  tropopause  showed  horizontal  uniformity  of  height  and 
backscatter  cross  section,  and  the  radiosonde  data  taken  near  the  radar  and 
near  the  common  volume  showed  similar  wind  and  temperature  structure  near 
the  tropopause,  the  signal  strength  on  the  forward  scatter  link  and  its  depen¬ 
dence  on  scattering  angle  behaved  in  accordance  with  the  prediction  of  turbulent 
scattering  theory  using  the  radar  data  as  an  input. 

The  radar  observations  have  shown  that  on  each  day  measurements  were 

made,  layers  were  detected  near  and  above  the  tropopause.  Turbulent  layers 

in  the  stratosphere  have  been  detected  at  heights  up  to  22  km.  These  layers 

provide  one  of  the  mechanisms  for  weak,  long-distance  troposcatter  propagation. 
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Monostatic  and  Bistatic  Scattering  from  Thin  Turbulent 
Layers  in  the  Atmosphere 

One  of  the  mechanisms  proposed  for  long  distance  transhorizon  propaga¬ 
tion  is  scattering  from  atmospheric  turbulence.  The  turbulent  scattering 
model  has  long  been  with  us.  The  usual  application  of  this  model  to  the  trans¬ 
horizon  scattering  problem  has  met  with  little  success.  The  dependence  of 
long  term  median  propagation  loss  measurements  on  frequency  and  on  scat¬ 
tering  angle  has  not  agreed  with  the  predictions  of  the  turbulent  scattering 
model.  Because  meteorological  data  point  to  the  existence  of  turbulent  regions 
in  the  atmosphere  and  because  of  a  requirement  for  a  technique  to  predict 
transhorizon  field  strengths  given  the  meteorological  data,  an  investigation  of 
scattering  by  atmospheric  turbulence  was  performed. 

Turbulent  Scattering:—  A  basic  reference  to  scattering  by  turbu¬ 
lence  is  the  work  by  Tatarski.  ^  In  this  work  the  bistatic  scattering  cross  sec¬ 
tion  per  unit  volume  is  related  to  the  spatial  power  spectral  density  of  the 
distribution  of  refractive  index  by 

/3i  (i?)  =  8Trk^  $  (k  -  km) 

/3|l  (t?)  =  /3j_ (t?)  cos2t? 

where  =  bistatic  scattering  cross  section  per  unit  volume  for 

polarization  perpendicular  to  the  plane  of  scattering. 
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j3||  (t?)  =  bistatic  scattering  cross  section  per  unit  volume  for 
polarization  parallel  to  the  plane  of  scattering, 
i?  =  scattering  angle. 

k  =  2tt/X  =  wave  number,  \  =  wavelength, 
k  =  vector  pointing  in  the  direction  of  propagation  with  a 
magnitude  equal  to  k. 

m  =  unit  vector  pointing  from  the  scatterer  to  the  receiver. 

$  =  spatial  power  spectrum  of  refractive  index  inhomogenities . 

n 

For  turbulent  refractive  index  fluctuations 

-  li/3 

9  (x)  =  $  (2k  sin  -J)  =  0.  033  C  Z  (2ksin  -?)  ;  .01  <— <  100  m 

n  n  2  n  2  x 

as  given  by  the  semi-empirical  similarity  theory  of  turbulence  for  the  inertial 

subrange.  The  classical  application  of  turbulence  scattering  theory  to  the 

2 

troposcatter  problem  has  been  to  assume  that  is  constant  over  the  entire 
common  volume  determined  by  the  intersection  of  the  two  antenna  patterns 
and  in  general  is  constant  over  all  space  or  decreases  exponentially  with 
height.  With  these  assumptions,  the  bistatic  scattering  cross  section  per  unit 
volume  can  be  extracted  from  measured  data  and  its  frequency  and  angular 
dependence  ascertained  through  the  parameter  x.  Figure  1  shows  the  depen¬ 
dence  of  the  turbulent  scattering  cross  section  upon  the  scale  size  i  =  2tt/x 
and  wavelength  for  three  values  of  the  scattering  angle  (referred  to  as  tropo¬ 
scatter  wavelength  where  the  scattering  angle  is  determined  by  the  geometry 
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of  the  troposcatter  link).  In  this  figure  are  shown  three  turbulent  scattering 

laws  each  for  a  different  region  of  i  .  The  scale  sizes  that  mark  the  extent 

of  a  region  are  not  as  well  defined  as  shown  in  the  figure  and  the  width  of  the 

different  regions  depends  upon  the  mean  properties  of  the  medium.  Figure  1 

shows  that  to  apply  turbulent  scatter  theory  for  the  inertial  subrange,  both 

the  frequency  and  geometry  have  to  be  properly  selected. 

A  second  problem  in  the  comparison  of  a  turbulent  scattering  theory  with 

measurements  arises  from  the  size  of  the  turbulent  region.  Radar  back- 

scatter  measurements  and  aircraft  measurements  show  that  the  turbulent 

scattering  volumes  typically  occur  as  thin  horizontally  stratified  layers  or 

patches.  These  layers  are  often  only  a  few  tens  of  meters  thick  and  tens  of 

kilometers  across.  An  example  of  a  radar  observation  of  a  thin  layer  is 

shown  in  Fig.  2.  From  atmospheric  model  computations  based  upon  the 

2  2 

methods  of  Vasil’ chenko  that  attempt  to  relate  the  C  value  to  the  mean 

n 

gradients  of  temperature,  humidity,  wind  speed  and  wind  shear,  it  has  been 

2 

shown  that  the  relative  magnitude  of  is  a  couple  of  orders  of  magnitude 
higher  in  the  layer  than  in  the  region  above  or  below  the  layer.  This  implies 
that  a  common  volume  averaged  bistatic  scattering  cross  section  per  unit 
volume  can  give  misleading  results  in  an  attempt  to  compare  measured  cross 
sections  with  theory.  It  also  implies  that  model  computations  of  transhorizon 
field  strengths  using  a  turbulent  scattering  model  should  use  a  volume  common 
only  to  the  layer  and  the  antenna  patterns.  Finally,  since  the  scattering  layer 
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is  thin,  the  layer  may  act  as  a  partially  reflecting  surface  for  vertical  scale 
sizes  larger  than  several  tens  of  meters  (frequencies  below  about  600  MHz 
on  a  190  km  scatter  link). 

Measurement  Program:—  A  scattering  model  based  upon  turbu¬ 
lent  scattering  theory  and  the  existence  of  isolated  thin  turbulent  layers  or 
patches  can  explain  away  many  of  the  objections  to  the  use  of  turbulent  volume 
scattering  as  a  transhorizon  scattering  mechanism.  The  questions  then  arise 
as  to  the  occurrence  of  these  thin  scattering  layers  and  a  demonstration  of 
transhorizon  scattering  as  caused  by  these  layers.  To  investigate  both  of 
these  questions  a  measurement  program  was  initiated  at  M.  I.  T.  Lincoln 
Laboratory  to  detect  high  altitude  layers  with  a  radar  and  correlate  their 
existence  with  the  results  of  long  distance  troposcatter  measurements.  The 
measurement  program  used  both  the  Millstone  Hill  L-band  radar  and  an 
X-band  troposcatter  link  operated  from  Wallops  Island,  Virginia  to  the 
Westford  Communications  Terminal.  The  common  volume  of  the  bistatic 
scattering  system  was  not  close  enough  to  the  radar  to  provide  a  direct  mea¬ 
surement  of  the  layered  structure  in  the  common  volume.  The  structure 
within  the  volume  was  inferred,  however,  from  the  radar  data,  radiosonde 
soundings,  and  the  angular  dependence  of  the  bistatic  scattering  data. 

The  characteristics  of  the  L-band  radar  when  used  for  weather  studies 

3 

are  mentioned  elsewhere.  The  radar  system  has  a  narrow  bandwidth  (0.  6°), 
a  relatively  large  pulse  length  (5  or  10  p  sec),  and  is  capable  of  detecting 
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turbulent  volumes  large  compared  with  the  resolution  volume  with  values 

higher  than  10  ^  m  (5  x  10  ^  cm  at  distances  of  100  km.  The 

radar  wavelength  is  0.  23  meters  and  for  turbulence  measurements,  it  selects 

a  scale  size  in  the  inertial  subrange  as  shown  in  Fig.  1. 

High  altitude  turbulent  scattering  as  measured  by  the  radar  is  generally 

weak  and  noise-like.  The  detection  of  turbulent  scattering  therefore  consists 

of  recognizing  a  spatial  pattern  of  noise-like  returns  embedded  in  the  receiver 

noise  when  displayed  on  a  range,  height  display.  Figure  2  shows  such  a 

display.  The  horizontal  lines  are  the  thin  turbulent  layers,  the  constant 

range  targets  are  from  ground  clutter,  and  the  curved  lines  are  aircraft 

2 

targets  detected  through  the  side  lobes  of  the  antenna  pattern.  The  value 

for  this  layer  is  obtained  from  post-test  analysis  of  50  pulse  incoherently 

averaged  received  signal  values  that  are  digitized  and  recorded  on  magnetic 

2  -16  -2/3 

tape.  For  this  particular  layer  =  5x10  m  at  a  point  10  km  high 

and  85  km  in  range.  The  true  value  may  be  higher  if  the  layer  is  thinner  than 
the  resolution  volume  of  the  radar. 

The  X-band  scatter  link  consisted  of  a  transmitter  located  at  the  JAFNA 
radar  facility  on  Wallops  Island  and  a  receiver  at  the  Westford  Communications 
Terminal  as  shown  in  Fig.  3.  The  transmitter  facility  consisted  of  a  1  kw  cw 
source  operated  at  7.  74  GHz  and  a  6-foot  linearly  polarized  antenna.  The  re¬ 
ceiver  facility  consisted  of  a  60-foot  circularly  polarized  antenna,  a  cooled 
parametric  rf  amplifier,  and  a  narrow  band  (500  cycle)  phase  lock  tracking 
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receiver.  The  minimum  detectable  signal  with  this  system  was  -160  dbm 
and  corresponds  to  a  troposcatter  path  loss  of  148  db  relative  to  free  space. 

The  great  circle  path  scatter  link  geometry  is  shown  in  Fig.  4.  The  path 
covered  a  great  circle  arc  length  of  628  km  and,  for  an  effective  four-third 
earth  radius  atmospheric  refraction  model,  the  height  to  the  lower  edge  of  the 
common  volumes  was  5.  5  km.  The  path  was  selected  so  only  the  regions  of 
the  atmosphere  above  6  km  were  illuminated  under  normal  propagation  condi¬ 
tions.  This  was  done  because  previous  radar  measurements  had  shown  that 
typically  one  to  three  isolated  scattering  layers  would  be  present  above  this 
height  and,  in  order  to  measure  the  weak  signal  scattered  by  these  layers, 
the  stronger  lower  level  layers  should  not  be  simultaneously  visible  to  both 
the  transmitter  and  receiver.  If  the  lower,  stronger  layers  were  visible,  the 
signals  scattered  from  these  layers  could  mask  those  of  interest  even  if  the 
strong  layers  were  in  the  side  lobe  patterns  of  both  antennas. 

The  forward  scatter  measurements  were  conducted  by  fixing  the  trans¬ 
mitter  antenna  pointing  direction  at  a  0.  8°  elevation  angle  and  along  the  great 
circle  azimuth.  The  receiving  antenna  was  scanned  in  azimuth  and  elevation 
so  as  to  include  the  region  of  the  layer  illuminated  by  the  transmitter.  The 
geometry  and  frequency  for  the  scatter  experiment  were  chosen  so  the  turbu¬ 
lent  scattering  process  is  described  by  the  above  equation  valid  in  the  inertial 

2 

subrange.  For  model  computations,  both  a  value  and  a  layer  thickness 
must  be  chosen.  Figure  5  shows  that  results  of  model  computations  for  layers 
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at  three  heights  and  for  elevation  scans  of  the  receiving  antenna  along  the 

great  circle  direction.  Although  the  scattering  angle  dependence  of  turbulence 

in  the  inertial  subrange  is  included  in  the  model  calculations,  the  shape  of  the 

curve  depends  primarily  on  the  transmitter  antenna  pattern  and  the  thin 

layered  geometry  of  the  scattering  volume.  Other  volume  scatterers  such  as 

the  ice  crystals  of  a  cirrus  cloud  will  give  curve  shapes  close  to  these  because 

the  scattering  angle  changes  only  slightly.  From  measurements  corresponding 

to  the  model  computations  illustrated  in  Fig.  5  only  the  layer  height  may  be 

2 

determined.  A  value  of  times  the  layer  thickness  may  also  be  determined 
for  the  layer  but  positive  identification  of  the  layer  as  turbulent  scattering  in 
contrast  to  particulate  scattering  cannot  be  made  from  these  elevation  scans 
alone.  Figure  5  was  generated  without  including  the  attenuation  due  to  gaseous 
absorption.  For  this  frequency  the  attenuation  is  of  the  order  of  3  db  and  is 
nearly  constant  for  the  range  of  layer  heights  and  angles  used. 

Results  of  Measurements:—  Forward  scatter  measurements  were 
made  on  19  days  during  the  months  of  March  through  May.  On  two  of  these 
days,  rain  along  the  path  provided  excessive  attenuation  and  signals  were  barely 
detectable.  The  signal  levels  for  the  remaining  days  ranged,  on  average, 
from  -144  dbm  to  noise  for  the  elevation  and  azimuth  angles  used.  Signals  in 
excess  of  -110  dbm  were  often  observed  due  to  aircraft  in  the  common  volume 
but  these  data  values  were  discarded.  Elevation  scans  were  made  on  14  of  the 
days  and  on  13  of  these  identification  of  a  primary  scattering  layer  could  be 
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made. 


Sample  data  for  the  13th  and  14  of  March  are  shown  in  Fig.  6.  On  the 

2 

curves  are  superimposed  the  best  fit  layer  model  and  the  to  which  it 

2 

corresponds.  The  values  are  for  an  assumed  layer  thickness  of  100  m 
and  may  deviate  from  this  value  due  to  the  actual  layer  thickness.  Because 
the  form  of  the  received  power  vs.  elevation  angle  curve  changes  significantly 
with  layer  height,  a  layer  height  may  easily  be  assigned.  For  both  days  the 
layers  are  near  8  km  in  height  and  agree  with  the  position  of  the  tropopause 
as  dete-  nined  from  radiosonde  soundings  made  close  to  the  center  of  the 
common  volume  (JFK  station  on  Fig.  2). 

A  direct  comparison  of  radar  and  forward  scattering  measurements  is 
impossible  because  of  the  distance  from  the  radar  to  the  common  volume.  The 
radar  data  has  shown  however  that  the  turbulent  patches,  when  they  occur, 
usually  occur  at  the  same  height  over  periods  of  several  hours  and  distances 
of  hundreds  of  kilometers.  The  actual  life  history  of  a  single  turbulent  patch 
may  however  be  brief  being  of  the  order  of  magnitude  of  tens  of  minutes. 
Figures  7  and  8  show  radar  detected  layers  at  heights  of  4,  6,  10  and  12  km 
for  8  May  1968.  The  data  shows  that  the  layered  regions  exist  over  several 
hundred  kilometers  near  the  Millstone  Radar.  If  we  assume  this  layer  struc¬ 
ture  extends  another  hundred  kilometers  into  the  common  volume  of  the 

2  -16 

scatter  system,  a  layer  height  of  10  km  and  with  ~  4  x  10  would  be  ex¬ 
pected.  The  lower  layer  of  those  in  the  common  volume  is  expected  since,  at 
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the  lower  elevation  angles  and  for  equal  scatter  cross  sections,  the  lower 
layer  produces  a  higher  signal  strength.  The  forward  scatter  measurements 
for  the  same  time  are  given  in  Fig.  9  and  clearly  show  that  the  primary  scat¬ 
tering  layer  is  at  10  km  and  is  in  agreement  with  our  expectation  from  the 
radar  data. 

Conclusions:—  A  thin-layered  volume  scattering  model  is  pro¬ 
posed  as  a  mechanism  for  long  distance  troposcatter  propagation  at  micro- 
wave  frequencies.  Using  a  turbulent  scattering  model,  the  expected  elevation 
angle  dependence  of  received  signals  scattered  from  single  layers  near  the 
tropopause  was  computed.  The  forward  scatter  measurements  on  13  of  14 
days  are  consistant  with  the  model.  For  days  with  radar  data  and  on  which 
radiosonde  data  show  similar  soundings  both  at  the  radar  and  at  the  common 

volume,  the  height  of  troposcatter  detected  layers  agrees  with  the  height  of 

2 

the  radar  detected  layers.  The  C  value  detected  by  both  the  radar  and  the 

n 

scatter  system  are  of  the  same  order  of  magnitude  when  compared.  The  dis¬ 
crepancies  between  these  estimates  arise  because  of  the  thinness  and  patchi¬ 
ness  of  turbulent  layers. 
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ft  (Relative  Units) 


Radar  Wavelength  (m) 

h — Wallops  -  West  ford  Link 
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Fig.  1.  Relative  scattering  cross  section 
per  unit  volume  vs.  wavelength  and  scale. 
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Fig,  2.  Millstone  Radar  measurements. 
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Fig.  3.  Coastal  map  showing  Wallops  Island  to  Westford  path. 
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Fig.  4.  Wallops  Island  -  Westford  cross  section. 
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Fig.  5.  Expected  variation  of  received  power  with  receiver 
antenna  elevation  angle. 
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Fig.  6.  Received  signal  vs.  elevation  angle  of  receiver  antenna. 
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Fig.  7.  Millstone  Radar  measurements  -  8  May  1968. 
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Fig.  8.  Millstone  Radar  measurements  -  8  May  1968. 
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Fig.  9.  Received  signal  strength  vs.  elevation  angle 
of  receiving  antenna. 
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